Probing the nanohydrodynamics at liquid-solid interfaces using thermal motion 
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We report on a new method to characterize nano-hydrodynamic properties at the liquid/solid 
interface relying solely on the measurement of the thermal motion of confined colloids. Using 
Fluorescence Correlation Spectroscopy (FCS) to probe the diffusion of the colloidal tracers, this 
optical technique -equivalent in spirit to the microrheology technique used for bulk properties- 
is able to achieve nanometric resolution on the slip length measurement. It confirms the no-slip 
boundary condition on wetting surfaces and shows a partial slip b — 18 ± 5 nm on non- wetting ones. 
Moreover, in the absence of external forcing, we do not find any evidence for large nano-bubble 
promoted slippage on moderately rough non- wetting surfaces. 

PACS numbers: 68.08.-p, 68.15.+e,47.15.Gf 



o 



i 

C 

O 

o 



> 

in 
o 

o 
in 
o 



T3 
C 

o 
o 



Over the recent years the pursuit of scale reduction 
inherent to nanotechnologies has been extended to the 
fluidic domain and liquid flow manipulation, with the 
important development of micro- and nanofluidics 0. 
However, reducing the scale of any system leads invari- 
ably to an enhancement of the influence of surface proper- 
ties with respect to the bulk ones: given the scale reduc- 
tion, most phenomena take place at the boundaries and 
a fundamental understanding of how surface properties 
might affect the overall flow properties has become cru- 
cial to design and optimize operational devices Q, 0, ■ 

Classically, one accounts for the influence of these in- 
terfaces through effective boundary conditions (BC) in 
the description of macroscopic hydrodynamics, the most 
common of those BC being the no-slip assumption (see 
Ref . |j| for an exhaustive review of the litterature) . How- 
ever, the possible deviation from this classical hypothe- 
sis, resulting in liquid slippage at the solid surface, has 
recenlty become a central issue, with immediate perspec- 
tives in the micro- and nanofluidics domains or in the 
electrokinetic context [3|- Slippage is usually accounted 
for by replacing the no-slip BC for the tangential veloc- 
ity v t by a partial slip BC, in the form bd z v t = Vt (with 
z perpendicular to the planar surface). This generalized 
BC introduces an extrapolation length b, usually denoted 
as the slip length 0,0. 

Probing the interfacial dynamics of liquids close to 
solid substrates has accordingly opened new experi- 
mental challenges in the recent years: specifically de- 
vised experimental tools capable to investigate the nano- 
hydrodynamics close to the solid substrate have been de- 
veloped, allowing to address the existence and conditions 
for slippage. Most recently, two main routes have been 
followed [4| : (i) dissipation methods on one hand on 
the basis of Surface Force Apparatus -SFA- and Atomic 
Force Microscopy -AFM- measurements; (ii) flow char- 
acterization close to surfaces on the other hand, using 
optical methods, like FRAP in evanescent waves geome- 
try, or microPIV velocimetry. These various approaches 



suggest an overall link between the hydrodynamic slip- 
page and the wettability of the solid substrate, in agree- 
ment with theoretical and numerical results 0, |fj . They 
fail however to provide a unified description of the slip- 
page phenomenon, reporting for instance slip lengths b 
on seemingly identical smooth surfaces that span several 
orders of magnitude, from tens of nanometers to microm- 
eters 0. The presence of nanobubbles at the surface, 
whose existence might be promoted by the driving flow, 
has been put forward to rationalize these results 0,0 E} • 
As demonstrated theoretically the existence of 

such nanobubbles would strongly enhance the measured 
(apparent) slip length, and could also provide an expla- 
nation for shear rate dependent effects reported together 
with large slip In order to clarify the experimen- 

tal picture, new non-intrusive approaches to probe the 
hydrodynamic surface properties are needed. 

In this letter, we demonstrate a completely differ- 
ent experimental route, allowing to explore the nano- 
hydrodynamics of liquids close to surfaces in the absence 
of any external forcing, thereby avoiding problems inher- 
ent, in all previous measurements, to the imposed flow. 
Rather than measuring the interfacial dissipation or the 
forced surface flow, we take advantage of the informa- 
tion already included in response to thermal fluctuations 
|l.'"{|| to extract the interfacial dynamics. This technique, 
which is analogous -for surfaces- to the passive microrhe- 
ology technique for bulk characterization [14j . proves to 
be extremely sensitive as it allows us to reach an un- 
precedented resolution for an optical technique, namely 
a few nanometers on the slip length measurement. To- 
gether with confirming the emerging picture for smooth 
surfaces of a moderate slippage (b = 18 nm) present only 
in non-wetting situation, it provides a truly non-intrusive 
technique to explore the possible role of gas pockets 
(nanobubbles) on slip properties. As a first step, it shows 
that in the absence of forcing, nano-bubble promoted gi- 
ant slip is absent on moderately rough non-wetting sur- 
faces. 
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FIG. 1: a) Scheme of the experimental cell with FCS detec- 
tion. Fluorescent colloids suspended in water are confined 
between a plane and a spherical lens. FCS : (MD) Dichroic 
mirror, (Obj.) microscope objective, (PH) pinhole, (APD) 
avalanche photodiode. b) Typical auto-correlation curve to- 
gether with theoretical fit according to eq. (0. 



We first describe the general principle of our approach. 
The diffusion dynamics of colloidal tracers is measured in 
a confined geometry between two solid surfaces of inter- 
est, using a home-built Fluorescence Correlation Spec- 
troscopy (FCS) device (fig. P). Tracer dynamics is af- 
fected by confinement and this dependence reflects the 
hydrodynamic boundary conditions which apply on both 
solid substrates 0, flil Il7| . Results for different nature 
of the solid substrates (in particular varying wettability 
and surface roughness) lead to measurable differences in 
the diffusion coefficient, allowing us to deduce the corre- 
sponding surface slippage. 

Let us now enter into the details of the experimental 
setup. Colloidal tracers (polystyrene [Molecular Probes] 
or silica [Kisker] fluorescent beads with typical diame- 
ters 2a ~ 200 nm and concentration c = lbead//im 3 ) in 
aqueous solution (10~ 5 M NaOH; 8 10 _3 M KC1) are con- 
fined between two solid silica surfaces made from a BK7 
spherical lens in contact with a Pyrex plane (fig.[TJ. The 
thermal diffusion dynamics of the colloids is measured in 
this confined geometry with a FCS device (fig-CJ. With 
this technique, beads fluorescence is excited with an ar- 
gon laser (Lasos) focused with a long working distance 
microscope objective (Leica x40, NA 0.8). Fluorescence 
is then collected through the same objective and sent to a 
detector (APD, Perkin- Elmer) connected to a correlator 
(Correlator.com©) via a dichroic mirror and a bandpass 
filter (Omega). While in bulk measurements, a confo- 
cal pinhole is inserted in the detection pathway to get a 
spatially defined measurement volume v, here the axial 
limits are set in pratice by the two confining walls so that 
v = (7rw 2 )e, with w the beam waist radius and e the wall 
to wall distance. With a bead concentration such that 
the mean number of tracers in the detection volume is 
low (typically around 1 here) fluctuations of the collected 
intensity I(t) arise due to the motion of beads entering 
or leaving the measurement volume. The characteristic 
time scale for such fluctuations corresponds to the resi- 
dence time td = w 2 / D of a bead within volume v, where 



D is the bead self-diffusion coefficient. More quantita- 
tively, considering a gaussian radial intensity distribution 
for the illuminating laser beam, the fluorescence intensity 
auto-correlation function reads 
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from which the experimental average number of beads n 
and their residence time can be extracted, as shown in 
fig-ffl As the bead diffusion coeficient D depends on the 
location z within the gap, we eventually obtain an av- 
eraged value over the all accessible positions within the 
thickness e. To ensure good statistics, values for each 
confinement e/2a were accumulated from over 100000 
events (associated with one bead passage) splitted be- 
tween many different auto-correlation runs and from at 
least 5 different experimental cells. Runs belonging to 
different cells may slighlty differ in absolute transit time 
and average number of beads due to small day to day 
variations in temperature, concentration or focus loca- 
tion. We accordingly normalized data from each cell by 
the reference point for (almost) unconfined beads, located 
at e/2a = 16.8, for which excellent statistics was already 
achieved in the time frame of a single cell experiment. 

For any location of the sphere-plane geometry where 
we measured n and 7~d, the wall to wall distance e was 
simultaneously measured using the interference pattern 
generated by the two confining surfaces (Newton rings). 
Due to the chosen bead size (2a ~ 200 nm), we were able 
to restrict our data point collection to the dark fringes 
(setting Ae = X/2n w = 183 nm; with A the laser wave- 
length and n w the optical index of water), therefore opti- 
mizing the signal to noise ratio. Owing to the large radii 
of curvature of the lenses used (from 250 to 500 mm), e is 
constant to better than 0.4 % over the measurement spot 
size 2w ~ 1 fim. 

The first element we focused on is the evolution of 
the mean bead number with the confinement defined as 
e/2a. This evolution is shown in figure [21 (inset) where 
we recover a linear behavior for the averaged number of 
beads in the volume v as expected: n — (c)(-KW 2 )e. This 
measurement provides an important check that no de- 
pletion or adsorption due to bead-surface interaction is 
detectable in our system [27}. Note however that the 
linear evolution does not extrapolate to e/2a = 1 as ex- 
pected for the sole excluded volume effect but to a slightly 
higher value (1.2 in the inset of fig- El> - This is associated 
with the additionnal "excluded volume" resulting from 
electrostatic repulsions between the wall and the beads. 
For silica walls, under the present conditions (electrolytes 
concentration), the excluded region is D C xd. — 20 nm 
thick (for both polystyren and silica beads), while it is 
found to be thinner -D G xci. — 13 nm for silica beads close 
to silanized walls. Such a distance agrees with the es- 
timate obtained from the balance between the thermal 
energy and the electrostatic wall-bead repulsion energy 
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FIG. 2: Evolution of the diffusion time with the confine- 
ment e/2a between hydrophilic walls for various colloidal trac- 
ers: (•) Silica beads (2a = 218 nm); (o) Polystyren beads 
(2a — 210 nm); ( — ) Theoretical prediction using Femlab 
and assuming no-slip at both walls. Inset: Evolution of the 
bead number with the confinement, (o) Polystyren beads 
(2a = 210 nm); ( ) linear regression. 



(calculated using the experimental Debye length and typ- 
ical zeta potentials for silica). 

We now come to the measurement of the residence time 
-or reciprocally to the diffusion coefficient- of the beads 
as a function of the confinement e/2a. When the confin- 
ing walls are hydrophilic, we expect a usual no-slip BC 
to apply at their surface H,0|. In such a situation, the 
bead mobility should be strongly reduced by the walls 
proximity [19j as was already verified experimentally by 
a few groups [20L l2lj . This wetting configuration was 
therefore the starting point from which we have explored 
the influence of surface properties on the hydrodynamics. 
Our results obtained in such situation (aqueous solution 
with silica walls) are presented in figure [21 together with 
theoretical predictions assuming no-slip boundary condi- 
tion on both walls. For low confinements e/2a 3> 1, a 
very sound, approximate solution for the mobility can be 
constructed on the basis of the Faxen approximate solu- 
tion for one wall 0, by adding the independent contri- 
bution of each wall to obtain the theoretical bead mo- 
bility |16(. However to avoid theoretical approximations, 
we conducted numerical resolution of the Stokes equa- 
tion for a (non-slipping) sphere moving parallel to the 
two confining solid walls (each characterized by a [pos- 
sibly different] slip length b). A finite element method 
was implemented using Femlab© (see in [l6j for details). 
This theoretical prediction for the mobility was then av- 
eraged over accessible z (from a + D eKC \, to e — a — -D eX cl.)- 
As is evidenced in figure [3 the agreement between ex- 
periments and theory with no-slip BC on the walls is 
excellent up to the strongest confinements. It provides 
another evidence supporting the fact that no-slip applies 
on our smooth wettable substrates 0, 0] . In addition, 
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FIG. 3: Evolution of the diffusion time of silica beads 
(2a = 218 nm) with confinement a hydrophobilic silica lens 
and a planar surface with varying surface property: (•) 
smooth (< lnm peak-to-peak roughness), hydrophilic sur- 
face; (>) smooth hydrophobic wall (OTS coated silica); 
(□) rough (40 nm peak-to-peak roughness) hydrophobic wall 
(OTS coated silica). Theoretical predictions are obtained us- 
ing Femlab, assuming no-slip on the lens, and different BC on 

the plane: ( — ) no-slip; ( ) partial slip b = 18 nm; (— ■ — ) 

negative slip length b — —40 nm (see text for details). As a 
guide, the dotted line (■ ■ •) shows the theoretical prediction 
for a b — 100 nm slip length. 



the fact that beads with very different chemistry display 
the very same behavior is another indication that our sys- 
tem is free from specific surface-bead interactions, beside 
the measured electrostatic repulsion leading to -D exc i.- 

We now consider the influence of surface properties 
on the hydrodynamic BC. For that purpose the silica 
plane was covalently coated with hydrocarbon chains 
using octadecyltrichlorosilan (OTS). The resulting resi- 
dence time of silica beads is plotted for different confine- 
ment against the previous results for wetting (plain sil- 
ica) plane. When beads are confined enough (e/2a < 5), 
residence times appear to be systematically shorter close 
to hydrophobic surfaces than to hydrophilic ones. This 
effect is more easily captured when normalizing this 
evolution by the theoretical behavior in the absence of 
slip: -D CX p.(e/2a)/_D n o-siip(e/2a). A failure of the no-slip 
boundary condition should result in a departure of this 
ratio from 1 . This is indeed what is observed in the inset 
of figure |3] where the systematic trend described above is 
best evidenced. Moreover it is possible to render this de- 
parture accurately by introducing a finite slip length b in 
the theoretical calculations [l7| • The fitted behavior 
agrees remarkably well with the experimental data pro- 
viding a slip length on smooth OTS coated silica planes 
of b = 18 ± 5 nm. This value is in perfect agreement 
with published results on the exact same surfaces using 
a dynamic-SFA It demonstrates the ability of our 

aproach to characterize the nanohydrodynamics close to 
interfaces at zero shear stress together with its unreached 
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sensitivity for an optical method. 

This sensitivity was confirmed in a subsequent set 
of experiments, investigating the role of surface rough- 
ness on hydrodynamical properties of interfaces. Our 
method, relying only on thermal fluctuations, is par- 
ticularly adapted for such situations as it avoids possi- 
ble surface modifications (creation and/or coupling with 
nanobubbles) present in all previously reported methods 
where external forcing is required. Before the hydropho- 
bic OTS coating, the silica planes were first treated 
30min. with a piranha solution (1 vol. H2O2, 2 vol. 
H2SO4) leading to rough topography 22] with a peak-to- 
peak height (checked by AFM on 5 x 5 fim 2 scans) around 
40 nm, as compared to below 1 nm for untreated sur- 
faces. Even on non-wetting surfaces, roughness kills the 
slippage, leading to reduced beads mobility as compared 
with smooth non- wetting surfaces. Actually the mobility 
is found to be below the one for a no-slip BC. This results 
comes from the fact that the hydrodynamic boundary is 
shifted toward the top of the roughness while the inter- 
ference gap measurements provides only the average wall 
position. Therefore the measured gap e (against which 
data are plotted) is larger than the hydrodynamic gap 
thus resulting in a translation of the experimental data 
that might be interpreted as a negative slip length of 
— 40 ± 20 nm ■ Again this demonstrates the sensitivity 
of our technique, showing differences for surface proper- 
ties affecting the nanorheology on scales as low as tens of 
nanometers. An interesting conclusion deduced from our 
data is that in the absence of external forcing, only the 
"negative" effect of roughness on slippage is evidenced. 
No huge enhancement (reaching micrometric slip length 
p^ l possibly promoted by trapped gaz is observed. 

We have therefore shown that the measure of the 
thermal motion of colloidal tracers provides an alterna- 
tive and very sensitive method -achieving nanometric 
resolution- to address the nano-hydrodynamics of sim- 
ple liquids close to surfaces. This method is able to 
give information on interfacial dissipation without any 
external forcing, by exploiting the intimate fluctuation- 
dissipation link. Working in the strict zero-shear rate 
limit, avoids shear induced alteration of the surfaces 
(such as possible nucleation of nanobubbles). In the 
present experiments we find no-slip on wetting Pyrex sur- 
face, and a finite slip length of 18 nm for water on hy- 
drophobic OTS surfaces. These results are in full agree- 
ment with recent findings using a dynamic-SFA |18| and 
with numerical results [6j . Moreover a nanometric rough- 
ness has been shown to kill slippage even on a hydropho- 
bic surface. This absence of large slip efffect at zero shear 
rate would confirm that the - still debated - presence of 
nanobubbles might originate in the flow itself, as was al- 
ready conjectured j3, [llj, |24j . To answer this question, it 
would accordingly be interesting to use our equilibrium 
approach on previously sheared surfaces. Alternatively 
other types of surfaces might be considered, in particular 



super-hydrophobic surfaces for which very large slippage 
is expected [2(|. Work along these lines is in 

progress. 
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